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Abstract 

After the discovery of the Higgs boson at the LHC, it is natural to start the research program 
on the precision study of the Higgs-boson couplings to various standard model (SM) particles. We 
provide a generic framework for the deviations of the couplings from their SM values by introducing 
a number of parameters. We show that a large number of models beyond the SM can be covered, 
including two-Higgs-doublet models, supersymmetric models, little-Higgs models, extended Higgs 
sectors with singlets, and fourth generation models. We perform global fits to the most updated 
data from CMS, ATLAS, and Tevatron under various initial conditions of the parameter set. In 
particular, we have made explicit comparisons between the fitting results before and after the 
Moriond 2013 meetings. Highlights of the results include: (i) the nonstandard decay branching 
ratio of the Higgs boson is less than 22%; (ii) the most efficient way to achieve the best fit for the 
data before the Moriond update is to introduce additional particle contributions to the triangular- 
loop functions of H'j^ and Hgg vertices; (iii) the la allowed range of the relative coupling of HVV 
is 1.01 io!i4' which means that the electroweak-symmetry breaking contribution from the observed 
Higgs boson leaves only a small room for other Higgs bosons; (iv) the current data do not rule out 
pseudoscalar couplings nor pseudoscalar contributions to the H^^ and Hgg vertices; and (v) the 
SM Higgs boson provides the best fit to all the current Higgs data. 



I. INTRODUCTION 



It is of very high expectation that the observed particle at the Large Hadron CoUider 
(LHC) [1112] is the long-sought Higgs boson of the standard model (SM), which was proposed 
in 1960s P]. At the end of 2011, both the ATLAS and CMS ^ experiments at the LHC 
have seen some excess of events of a possible Higgs boson candidate in the decay modes of 
H ^ 'j'j, H WW* i~^ui~u, and H -> ZZ* U channels. Finally, the discovery 
was announced in July 2012 by ATLAS [I] and CMS [2]. The channels WW and ZZ are 
consistent with the predictions of the SM Higgs boson, while the 77 rate is somewhat higher 
than expectation. Some evidence is seen in the hh mode at the Tevatron [5], but the mass 
range is quite wide. On the other hand, the r+r^ mode appears to be suppressed before 
the very recent update, although the data contain large uncertainties. 

A previous update was presented at the Hadron Collider Physics Symposium 2012 P, [7] 
and in a series of experimental notes [5HT5] at the end of 2012. At that time, the r+r^ data 
began to appear, but still too early to say something concrete. The diphoton production 
rate was somewhat higher than the SM prediction by a factor of 1.4 — 1.8. Nevertheless, the 
deviations are only 1 — 2 a. A large number of models have been put forward to account 
for the observed particle at 125 GeV, including the SM, supersymmetric models such as the 
minimal supersymmetric standard model [15], the next-to-minimal supersymmetric standard 
model [T7], the f/(l)-extended minimal supersymmetric standard model [IB], fermiophobic 
Higgs boson [19j, two-Higgs-doublet models (2HDM) of various types [20j, Randall-Sundrum 
radion [2T] , inert-Higgs doublet model [22] , etc (a summary of various models can be found in 
Ref. [23].) They all can explain the enhanced diphoton rate with some choices of parameter 
space. Yet, more data are needed in order to firmly establish the excess in the diphoton 
channel. The most recent update was during the Moriond 2013 meetings [24] . The updated 
data can be found in a number of conference note from the ATLAS [251 |2S] and CMS [27H5U] . 

Currently, a number of decay and production channels are available. On the production 
side, there are gluon fusion (ggF), vector-boson fusion (VBF), and associated production 
with & V = W/Z boson (VH) and top quarks (ttH); while the decay channels include 77, 
ZZ* — 7- Ai, WW* — > i^ui"!?, bb, and r+r^. One can extract useful information on the 
size of the Higgs boson couplings from the available data. However, in order to do that 
the dependence of various production and decay modes on the Higgs couplings has to be 
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taken into account correctly. For example, the ggF depends on the Higgs couplings to a 
pair of top (Htt) and bottom (Hbb) quarks, as well as possible existence of exotic colored 
particles running in the loop; while the VBF and VH depend only on the Higgs coupling 
to a pair of vector bosons {HVV). Also, the decays into WW* and ZZ* simply depend on 
HVV, and the decays into bb and r+r~ depend on Hbb and Htt respectively; but the decay 
into 77 involves all of the above couplings and perhaps new electrically-charged particles 
in the triangular loop. A global analysis of all the Higgs couplings using all the available 
data would be extremely useful to identify the observed Higgs boson. Once we disentangle 
each of the Higgs couplings from the global data set, we can use the result to compare with 
models. This approach is in contrast to those top-down approaches, which usually start 
with a model, calculate the signal strengths, and then find the allowed parameter space to 
fit to the data. 

Indeed, the Higgs precision era just begins. There have been a number of works in the 
past few months going in this direction, in a more or less model-independent framework 
[3THSI], in 2HDM framework |H2HSS], and in supersymmetry [^0145^ . Also, there are studies 
toward the determination of the spin-parity nature of the Higgs boson (see Ref. [55] for more 
references in literature) that cannot be obtained from the signal strengths. 

About a couple of weeks after we posted the first version of our paper to arXiv, both 
ATLAS and CMS Collaborations have updated their Higgs data during the Moriond 2013 
meetings [21]. They have released a series of conference notes [261430] on the new data. In 
particular, the most striking is the change of the diphoton data by the CMS [27]. Because 
of this change the overall fits also change dramatically. In the following, we will show the 
results before and after the Moriond 2013 meetings. 

The characteristic features of our analysis are summarized as follows. 

1. We allow the Yukawa couplings to the charged-lepton (Htt) and the down-quark 
(Hbb) sectors to vary independently. This can be realized in some versions of the 
2HDM. This has also been adopted in a few previous works. 

2. We allow an independent deviation in the total decay width of the Higgs boson, in 
addition to the parameters of the Yukawa and HVV couplings. This can be realized 
in nonstandard decays of the Higgs boson, e.g., h — > XiXi in SUSY models, H ^ aa 
where a stands for some lighter Higgs bosons in the model. This has been included in 
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some previous works. An interesting result is obtained because of this improvement. 
The nonstandard decay width is constrained to be less than 1.2 MeV at 95% CL, 
which accounts for a branching ratio of about 22%. 

3. In the loop vertices of if 77 and Hgg, we allow new parameters AS"^ and AS*^ respec- 
tively, which can most conveniently account for the effects of new particles contributing 
to the loops. 

4. Preliminarily, the pseudoscalar interpretation of the observed Higgs boson is disfa- 
vored. However, it may not apply to the case of CP-mixed state carrying both scalar- 
type and pseudoscalar-type couplings simultaneously. We perform a whole new anal- 
ysis including both scalar-type and pseudoscalar-type Yukawa couplings, and both 
radiatively-induced scalar and pseudoscalar H'^'^ and Hgg vertices. We show that 
pseudoscalar couplings are actually not ruled out based on signal strengths only, by 
giving equivalently good fits compared to the CP-conserving case. 

The organization of the paper is as follows. In the next section, we describe the interac- 
tions of the Higgs boson, including deviations in the Yukawa couplings and deviations in the 
loop functions of -^77, Hgg, and HZ'-y vertices, as well as the notation used in the analysis. 
In Sec. HI, we list the Higgs boson data both before and after Moriond 2013 meetings that 
we use in this analysis. We present the results of various fits in Sec. IV, and the readers can 
see how the fits change when the Higgs data are changed. In Sec. V, we present the results 
using both scalar-type and pseudoscalar-type couplings. We conclude in Sec. VI. 

II. FORMALISM 
A. Higgs Couplings 

We follow the conventions and notations of CPsuperH [USHSH] for the Higgs couplings to 
the SM particles assuming the Higgs boson is a generally CP-mixed state without carrying 
any definite CP-parity. 

• Higgs couplings to fermions: 

^Hff = - E E ^/ + </775) / . (1) 

f=u,d,l ^^^^W i=i 
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For the SM couphngs, Q^jj = 1 and g^ff = 0- 
Higgs couphngs to the massive vector bosons: 
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Chvv = 9 Mw [9„^^W+W-^ + g„,, ^ Z,Z^^ H . (2) 

For the SM couphngs, we have g^^y^, = g^zz = 9hvv ~ respecting the custodial 
symmetry. 

Higgs couplings to two photons: The amplitude for the decay process — )■ 77 can be 
written as 

where ki^2 are the momenta of the two photons and ei,2 the wave vectors of the cor- 
responding photons, e^_|_ = — 2A;f(A;2 ■ e2j_ = £3 — 2/^2(^1 ■ £2)/^^ 
{eie2kik2) = e^yp^ 6162^1^2 • The decay rate of if — )■ 77 is proportional to |S"^p + |P'>'p. 
Including some additional loop contributions from new particles, the scalar and pseu- 
doscalar form factors, retaining only the dominant loop contributions from the third- 
generation fermions and W^, are given by[^ 

S\Mh) = 2 NcQ]g!jjfF,f{Tf)~g„^^F,{Tw) + /^S\ 

f=b,t,T 

P\Mh) = 2 Y: NcQ}g^fjF,f{Tf) + AP\ (4) 

f=b,t,T 

where = M^/4m^, Nc = 3 for quarks and Nc = 1 for taus, respectively The 
additional contributions AS"^ and AP^ are assumed to be real in our work, as there 
are unlikely any new charged particles lighter than Mh/^,. 

Taking Mh = 125.5 GeV, we find that 

5^ ~ -8.35 gnww + l-7Q9kt + (-0-015 + 0.017 z) gf^-,, 

+ (-0.024 + 0.021 i) gl^^ + (-0.007 + 0.005 i) g'^^^ + AS^ 

P'< ~ 2.18g'Hu + (-0.018 + 0.018 g^^u 

+ (-0.025 + 0.022 i) g^^^ + (-0.007 + 0.005 i) g^^^ + AP^ (5) 

giving S^^ = -6.64 + 0.043 i and Pg\^ = 0. 



^ For the loop functions of Fsf^pf^i{T), we refer to, for example, Ref. 
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• Higgs couplings to two gluons: Similar to — > 77, the amplitude for the decay pro- 
cess H ^ gg can be written as 

where a and b {a,b — 1 to 8) are indices of the eight SU{3) generators in the adjoint 
representation. The decay rate oi H ^ gg is proportional to j-S^p + |P^p. Again, 
including some additional loop contributions from new particles, the scalar and pseu- 
doscalar form factors are given by 

S'^iMn)^ Y.g'HfyF,f{Tf) + AS\ 

f=b,t 
f=b,t 

The additional contributions A^*^ and AP^ are assumed to be real again. 
Taking Mh = 125.5 GeV, we find that 

-5^ ~ 0.688 ^If, + (-0.037 + 0.050 i) gfj-,, + AS^ , 

~ 1.047^:^,-, + (-0.042 + 0.050Z) g^-,, + AP^ , (8) 

giving Sly^ = 0.651 + 0.050 i and PgM = 0. 

• Higgs couplings to Z and 7: The amplitude for the decay process H Z{ki, ei) ■y{k2, £2) 
can be written as 

Mz^H = - ^ {S^'{Mh) [h ■ k2 el ■e*-k,-e;k2- €{] - P^\Mh) {ele^kik^)} 

(9) 

where ki^2 are the momenta of the Z boson and the photon (we note that 2ki • k2 — 
Mfj — Mz), ei,2 are their polarization vectors. The scalar and pseudoscalar form factors 
are given by 

S'\Mh) = 2 E QfNlm'f ~ ^"^^^ gf,j^ + M| cot Owg^^^F^ + AS'^ , 

f=t,b,r ^WCW 

P'\Mh) = 2 E QM^) ~ ^'^^^ 9ln Pf + ^P"' ■ (10) 
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The additional contributions AS^'^ and AP^'^ are assumed to be real. The loop 
functions are El 



fP = Co{m})+AC2{m}), 
Ff = Co{m]) , 



-2c2 ) + 2(1 -6c2 



w) 



C2(M^) + 4(1 - 4c^)Co(M^) . (11) 



Taking Mh = 125.5 GeV, we find 



S^^ ~ -11.966 gnww + 0.615 gfj-^t + (-0.008 + 0.004i) g'^-^^ 
+ (-0.0004 + 0.0002 i) g%^^ + AS^^ , 

P^^ ~ Qm?>glu + (-0.009 + 0.0040 g^-,, 
+ (-0.0004 + 0.0002 i) g%^^ + AP^^ , 



(12) 



giving S^^ = -11.358 + 0.004 i and PglJ = 0- 

In passing, we recall that the Z-boson couplings to the quarks and leptons are given 
by 



(13) 



where gz = e/{swCw), v^ff = li/^ — QfS^ and a^// = -^1/2 with Jg = 1/2 and 



Jg^'' = -1/2. 



Finally, we define the ratios of the effective Higgs couplings to gg, 77, and Z7 relative to 
the SM ones as follows: 



C„ = 



\S9\^ +\P9\^ 



9 — \ I Cff |2 ' ^7 

> PsmI 



|^7|2 _^ |p7|2 



7 |2 ' 
SM 



C^7 



(14) 



Note that the ratios of decay rates relative to the SM are given by |Cgp, IC7P, and |C, 
respectively. 



^7! ) 



For the functions of C{),2{m?), we refer to 
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B. Signal strengths 



The theoretical signal strength may be written as the product 

fi{V,V)^fi{V)il{V) (15) 

where V — ggF, VBF, VH, ttH denote the production mechanisms and V — 77, ZZ, WW, bb, tt 
the decay channels. More explicitly, we are taking 





\S%MH)\' + \P3{MH)f 




/i(VBF) 


HWW,HZZ ' 


j5(VH) 


HWW,HZZ ' 


/2(ttH) 





and 



with 



(16) 



Note that we introduce an arbitrary non-SM contribution AFtot to the total decay width. In- 
cidentally, Ttot{H) becomes the SM total decay width when Gnff — 9Hff ~ ^1 9hwwhzz — 
1, AS'^'^'^'^ = AP"''^'^"' — 0. 

The experimentally observed signal strength should be compared to the theoretical one 
summed over all production mechanisms: 

fx{Q,V)^ Y: CQ-pn{V,V) (19) 

7'=ggF,VBF,VH,ttH 

where Q denote the experimentally defined channel involved with the decay V and the 
decomposition coefficients Cq-p may depend on the relative Higgs production cross sections 
for a given Higgs-boson mass, experimental cuts, etc. 
The associated with an uncorrelated observable is 
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where a^^^{Q,I)) denotes the experimental error. For two correlated observables, we use 



— + — 







[^EXP(g^^2))] 





;i-p^) 

(21) 



-2p 

where p is the correlation coefficient. 
C. Parameters using in the fits 

Without loss of generality we use the following notation for the parameters in the fits 

— UHuu ■> '^d — yudd 1 ~ 'dull ' ~ dnvv J 

— yHuu 1 ^d — iJudd ' 



P 

9mi 



(22) 



Here we assume generation independence and also custodial symmetry between the W and 
Z bosons. Note that the tree-level pseudoscalar couplings to W and Z bosons are zero. The 
ffist and second generation fermion couplings to the Higgs boson is rather small, but if in 
the near future the H — )■ p^ ji~ can be measured, one may set independent parameters 
and Cr (for the present work we consider them to be the same.) 
In the fits, we further use 



A5^ A^'^; AP^ AP^ 



(23) 



which are real quantities assuming that any new particles running in the triangular loop are 
heavier than one half of the Higgs boson mass. Note that the quantities and S"^ are in 
general complex in both the SM and beyond the SM. For the most direct comparison with 
experimental data we use and Cg in the plots. 

Another important quantity is the additional contributions to the width of the Higgs 



boson, AFtot as in Eq. (18). The parameter AFtot takes into account the nonstandard 
decays of the Higgs boson, e.g., hi — )■ aiai in NMSSM, hi — )■ XiXi in SUSY, /i2 — hihi in 



The quantities AC^ and ACg defined in Ref. [4^ are in general complex. It would be more difficult to 
use them as fitting parameters. 
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other extended Higgs models. The current data still allow a small amount of Higgs invisible 
decay branching ratio. 

We first use the following scalar-type couplings: 



C^, C^, Cf, a, /^S\ AFtot 

to fit to the Higgs data in Sec. IV, where we focus on the SM-like Higgs boson. In Sec. V, 
where we also consider the possibility that the observed Higgs boson can allow some level 
of pseudoscalar-type couplings 

Cj, , AP\ AP^ , 

in addition to the scalar ones. 

D. Two-Higgs Doublet Models 

Two-Higgs-doublet models employ two Higgs doublets in the process of electroweak- 
symmetry breaking (EWSB). A discrete Z2 symmetry is usually imposed in order to avoid 
dangerous tree-level fiavor-changing neutral currents. The most studied are the type I and 
type II models. They can easily be covered by the framework presented in this paper. We 
illustrate using the model II. 

The Higgs sector consists of two Higgs doublets = (H^ and — {h^ Hf^ 
where the subscripts d denote the right-handed quark singlet fields that the Higgs doublets 
couple to. After EWSB, there are two CP-even, one CP-odd, and a pair of charged Higgs 
bosons. The parameters of the model in the CP-conserving case can be chosen as 

f^hi f^H, f^A, tanp = — , a 

Vd 

where a is the mixing angle between the two CP-even Higgs bosons. The couplings of the 
fighter CP-even Higgs boson h (assumed to be the observed boson) to the tau, bottom, top 
quarks, and W/Z boson relative to their corresponding SM values are given by 

T-T+ hi tt W+W-/ZZ 

h: — sin a/ cos ^ — sin a/ cos ^ cos a/ sin ^ sin(/3 — a) 

Therefore, we can equate these quantities with the definitions of C^,C^, Cf , Cy, given by 
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From the above relations, one may derive the following consistency relations which should 
hold in the type II model: 

COS /3 = ^ ^ , sm /3 = — ; 

On the other hand, the only additional particle that can run in the triangular loop of /i77 
is the charged Higgs boson. Also, there are no new particles other than the SM particles 
that the Higgs boson h can decay into. Therefore, the other quantities 

A^^ ^ 0, AS^ = 0, AFtot = . (26) 

Thus, the two-Higgs doublet models, in general, can be covered by our framework. 

In more complicated Higgs sectors, e.g., with additional singlets, there may be lighter 
Higgs bosons that the observed Higgs boson can decay into. In such a case, the additional 
decay modes will contribute to AFtot- 



E. Models with singlet Higgs bosons 

Simple extensions of the SM Higgs sector with one or more Higgs singlet fields are attrac- 
tive, because they can often provide a dark matter candidate once some kinds of discrete 
symmetries are imposed on the extra fields. Some variants can be found in Refs. [7DH75] . 
In the simplest version 173J, the Higgs sector consists of the usual SM Higgs doublet $ and 
a real singlet Higgs field x- They couple to each other via a renormalizable interaction 
PX^ A discrete Z2 symmetry is imposed on % — i- — x such that x cannot develop the 
vacuum expectation value (VEV) and becomes a dark matter candidate. After $ develops 
the VEV, X couples to the H via the interactions x^H and x^H^. Therefore, the Higgs 
boson, in addition to the couplings to the SM fermions, also couples to a pair of xs- The 
only modification in our framework is the total decay width, accommodated by AFtot- 



F. Supersymmetric models 



There are many varieties in supersymmetric models which contain at least two Higgs 
doublets. In the minimal supersymmetric extension of the SM (MSSM), there are three 
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neutral Higgs states and, in principle, any of them can be the candidate for the observed 
particle at 125 GeV. If the i-th {i = 1,2,3) Higgs state is assumed to be the observed 
particle, we have 

C! = ^, C! = C^ = ^ a = 0^,cos/3 + 0^,,sin/3; 
sm jD cos jD 

Cf = -^0... Cf = Cf = -!^^0„, (27) 

sm p COS [3 

where O^^i^^^i and Oai denote the CP-even and CP-odd components of the z-th Higgs state, 
respectively |^ and cos /3 and sin /3 are defined in the same way as in the type-II 2HDM. 
Including the threshold corrections to the third-generation Yukawa couplings, the above 
tree-level relations undergo some changes. Nevertheless, this case can be covered because 

S P S P 

we are treating ' and ' independently in our framework. 

Beyond the MSSM there could be more than three neutral Higgs states. In this case, it is 



straightforward to find similar relations as in Eq. (27). For example, in the next-to-minimal 
supersymmetric standard model (NMSSM) the neutral Higgs bosons have a 5 x 5 mixing 
matrix [71] . 

Any bosonic and fermionic contributions of SUSY particles to the 7/77 and Hgg vertices, 
including the charged Higgs-boson contribution, can be nicely accommodated by using the 
parameters AS*^'"^ and AP^'*^. Also, the parameter AFtot can take into account any Higgs 
decays into the non-SM particles. 



G. Little Higgs models 

Little Higgs models belong to a class of models in which the quadratic divergences to 
the Higgs boson are cancelled by a set of particles having the same spin statistics as the 
SM particles. For each SM particle there corresponds a little-Higgs (LH) partner with the 
coupling to the Higgs boson specifically designed in such a way that the quadratic divergence 
is cancelled. For example, the W boson has the LH partner Wh- A phenomenological inter- 
esting example is the littlest Higgs model [7S], the phenomenology of which was described 
in details in Ref. [75]. 

In general, the Yukawa couplings of the SM-like Higgs boson could be different from the 
SM, depending on the gauge structure of the LH model; so are the couplings to the W/Z 
^ For the precise definition of the orthogonal 3x3 Higgs-boson-mixing matrix O, we refer to Ref. [55] , 
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bosons. Other heavy LH partners, if they are electrically charged, can run in the triangular 
loop of if77 vertex, and if they carry color they will contribute to Hgg vertex. If there are 
other light Higgs bosons that the SM-like Higgs boson can decay into, they will increase the 
decay width Ftot- Thus, the LH models can be accommodated in the present framework by 



Recent analyses of little Higgs models with respect to current Higgs data can be found 



H. Fourth Generation model 

The sequential fourth generation model is a simple extension of the SM by adding an 
analogous repetition of a generation of fermions. The new charged leptons and quarks can 
run in the loop of 1/77 while the colored quarks run in the loop of Hgg. Thus, the fourth 
generation contributes to AS"^ and AS*^ only while all the Yukawa couplings are the same 
as the SM and AFtot = 0. 

III. HIGGS DATA 

Current Higgs data focus on a few decay channels of the Higgs boson: (i) h — > 77, (ii) 
h ZZ* i+ti+t, (iii) h WW* i^vtv, (iv) h bb, and (v) h t+t~. 
Within each decay mode both CMS and ATLAS have reported a number of channels, such 
as inclusive, vector-boson-fusion tagged, and/or VH tagged. All the available ATLAS, CMS, 
and Tevatron data in these five decay channels are shown in Tables [iffV] respectively. Both 
sets of data before and after the Moriond 2013 meetings are listed in the tables. We have 
used 22 data points in our analysis. The chi-square of all these 22 data points relative to 
the SM is about 17.5 and 18.94 for the data set before and after Moriond, respectively, and 
so the chi-square per degree of freedom (dof) is about 17.5/22 = 0.80 and 18.94/22 = 0.86, 
which means that the SM is a reasonably good fit to the Higgs data. The goodness of the 
fit, measured by the p-value, is about p = 0.74 and p = 0.65 for the data before and after 
Moriond, respectively. |^ 

^ Assuming the goodness-of-fit statistics follows a probability density function, the p-value for the hy- 
pothesis is given by [75] 



c: 



iS,P 



a, AS^, and AL- 



tot- 



in Refs. [771 [78]. 
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There are four production modes: gluon fusion (ggF), vector-boson fusion (VBF), asso- 
ciated production with a W/Z boson (VH), and associated production with a tt pair (ttH). 
The production cross sections for each production modes at the LHC could be found in 
Ref. [80]. For a/s = 7 TeV and Higgs-boson mass Mh = 126 GeV, the cross sections are: 

cT(ggF) = 15.080, a(VBF) = 1.211, a(VH) = 0.8653, (T(ttH) = 0.0843 pb, (28) 



where (t(VH) = a{WH) + a{ZH). For ^ = 8 TeV and Higgs-boson mass Mh = 126 GeV 
the cross sections are: 

(T(ggF) = 19.220, (t(VBF) = 1.568, a(VH) = 1.0625, a{tm) = 0.1271 pb. (29) 

Since the ATLAS and CMS from the above tables combined both ^/s = 7 TeV and 

a/s = 8 TeV, we take the luminosities for the 7 and 8 TeV data as weights to recalculate the 

cross section of each production mode. Take the VH mode of "untagged" channel of "CMS 

(5.1 fb"^ at 7 TeV + 19.6 fb^^ at 8 TeV)" (after Moriond 2013) from Table |l| as an example, 

the weighted cross section of the VH mode is 

,,,,,, 5.1 fb^i X 0.8653 pb + 19.6 fb"^ X 1.0625 pb , , , 

aiVRU^.Me, = 5.1 fb 1 + 19.6 fb'^ = ^'^^^ • 

For other production modes we get 

a (ggF) weighted = 18.365, a (VBF) ^etghted = 1-494, a (ttB.) weighted = 0.118 pb . (31) 



By using the weighted cross sections, we obtain the decomposition coefficients Cqp (19) for: 
(i) "untagged" channel of CMS in Table |l( 

(ii) "Inclusive" channels in Table [H] and Table |III| 

(Hi) "/i(VBF + VH,rr)" channel in Table [V| under the assumption that the ggF and ttH 
production modes do not contribute. 

For the decomposition coefficients for "0/1 jet" and "VBF tag" in Table [V] we take the 
results for the three search channels fiTh + X, cTh + X, and e/x + X, presented in Tables 1, 
2, and 3 of Ref.[l5]. 



where n is the degrees of freedom and 

^n/2-lg-z/2 



fiz;n) 



2"/2r(n/2) 
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For the CMS "VBF tagged" channel in Table|T|and the "0/1 jet" and "VBF tag" channels 



in Table III, we borrow the numbers found in Ref. |17] . 



The Tevatron decomposition coefficients in Tables [T] and III are from the ratios of the SM 



Higgs production cross sections. Note that we do not use the Tevatron rr data upon the 
large uncertainty recently reported in Ref. [81j . 

IV. CP CONSERVING FITS 

In the CP conserving fits, we have fixed 

C^Ai = ^P'^' = ' (32) 

while varying 

Cj,,„ a, A^^'^ AFtof (33) 
More precisely we have implemented the following 5 fits: 

A. SMfit. 

B. One-parameter fit varying AFtot with C^^^ = = I and AS^'^ = 0. 

C. Two/three-parameter fit varying AS^ and A^^ without /with AF^t taking C^^^ = 
a = 1. 

D. Four-parameter fit varying Cf^^ and with AS^'^ = and AFtot = 0. 

E. Six-parameter fit varying rf^, C^, and AS^''^ with AFtot = 0. 

As mentioned above we are going to show the fitting results with the data set before and 
after the Moriond 2013 meetings for each fit. Due to the change in the data, especially the 
diphoton data from the CMS, the fitting results change dramatically. In each of the fits, we 
first describe the fitting results with respect to the data after the Hadron Collider Physics 
Symposium 2012 but before the Moriond 2013; then the fitting results with respect to the 
data after the Moriond 2013. The best-fit values for the parameters of the above fits are all 
summarized in Table IVll 

The figures are shown for the fits with the most updated data after Moriond. Note also 
that when we show the 2D regions for two of the varying parameters in the figures, we 
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are marginalizing over the other parameters if more than two are allowed to vary. This also 
applies to the next section of CP violating fits. 

A. SM fit 

1. Before Moriond 

As we have mentioned the SM fit gives a /dof = 17.5/22 = 0.8, It gives a p-value of 
p = 0.74, which means the SM has a chance of 0.74 to be the true interpretation of the data. 
Contributions of chi-square from each data are shown in the second last column of Tables |I} 
|Vj The i7 — 7- 77 data from ATLAS, CMS, and Tevatron give the largest contribution to 
the chi-square, while the least contribution is from H — t- ZZ*. Specifically, in each decay 
channel the largest contribution is from: ATLAS figgH+uH and Tevatron {H — )■ 77); CMS 
Inclusive {H ZZ*); CMS VBF tag {H WW*); ATLAS VH tag {H hi); ATLAS 
fi{ggF, VBF+VH) {H — > rr). As we shall see soon that because the chi-square is dominated 
by the diphoton data the most efficient way to fit to the data is using AS"'' and AS^. 

2. After Moriond 

The SM fit gives a /dof = 18.94/22 = 0.86, which gives a p-value of p = 0.65. This 
value shows that the SM description of the data stays more or less the same as before the 
Moriond update. The diphoton data still dominate the total x^- The x"^ of each decay 
channel, shown in the last column of Tables |I]jV} is about the same as before. 

B. Vary only AFtot while keeping = C| = Cf = = 1 and AS^ = AS^ = 

1. Before Moriond 

We found that varying AF alone does not improve the chi-square. Numerically the chi- 
square per dof is 17.5/21 and the 95% allowed range for AFtot is 

-0.022+^:^^ MeV . 

The central value is consistent with zero and thus the 95% CL upper limit for AFtot is about 
1.4 MeV. Note that the total width of the SM Higgs boson with Mh = 125.5 GeV is about 
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4.1 — 4.2 MeV. Therefore, the 95% CL upper hmit for the nonstandard branching ratio of 
the Higgs boson is about 25%. The nonstandard decays include invisible decays, decays into 
other lighter Higgs bosons, or decays into other exotic particles. 

2. After Moriond 

The situation remains the same. The /dof = 18.89/21 and the 95% allowed range for 
AFtot is 

0.10+J:^^ MeV . 

The central value is consistent with zero and thus the 95% CL upper limit for AFtot is about 

1.2 MeV. Therefore, the 95% CL upper limit for the nonstandard branching ratio of the 
Higgs boson is about 22%. 

C. Vary AS^ and AS^ while keeping = = Cf = Cy = 1 

1. Before Moriond 

In this fit, only the parameters AS"'' and AS*^ vary, which is simply motivated by the fact 
that the most deviated Higgs data are the diphoton signal strengths while all the other data 
are more or less consistent with the SM values. It turns out that this is the most efficient 
way to fit to the data statistically (the most efficient here means that the is reduced the 
most with the best /dof .) The best fit values are 

A^^ = -2.731};}^, AS^ = -0.050^°:°^t xV^^o/ = 11.27/20 = 0.56 (34) 

There are two solutions of AS*^, which can be easily understood from the expression of 
5*^ in Eq. (|8]). Since the signal strengths depend on the absolute value of S*^, numerically, 
±0.6 ~ 0.65 + ASf, which gives AS^ ~ -0.05 or -1.25. They both give Cg ~ 0.92. Also, we 
were to extend the range of AS"'' further, there would have been a solution around AS'^ ^ 16, 
according to the expression of S"^ in Eq. ([s]). Nevertheless, it is unrealistic to find a model 
that can generate such a large AS"''. The quantities C-^ and Cg are very close to physical 
observables, and so their best values are 

- 1.41, Cg - 0.92 . 
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The chi-square per dof for this two-parameter fit is 0.56 compared to 0.8 of the SM, which 
shows some real improvement. 

If we further allow AFtot to vary simultaneously with AS^ and AS^, the best fit is 

= -2.93^J:^?, A^^ = 0.0063lE^:J^, AFtot = 0.791?;? MeV, xVdof = 10.83/19 = 0.57. 

It is clear that including AFtot in the fit does not improve the chi-square per dof. Since 
AFtot is still consistent with zero in this case, we will fix AFtot = in the later fits. 

2. After Moriond 

The most obvious difference between the data set before and after the Moriond 2013 can 
be seen here in this fit. Before the Moriond the is dominated by the diphoton data, in 
which both ATLAS and CMS showed 1.5 — 2a excesses, and so the dynamics of the fit will 
push to the direction to substantially reduce the of the diphoton data. However, with the 
new CMS diphoton data (0.78lo;26) the whole fit changes. The dynamics of the fit cannot 
force the parameters to go into one direction, because the ATLAS data is still about 1.5a 
larger than the SM while the CMS one is about la smaller. 

The best fit values for AS"'' and AS^ are 

AS'-' = -0.96l[J|^, A^^ = -0.043 ± 0.052, xVdof = 17.55/20 = 0.88 (35) 

We notice that the sizes of the errors are reduced reflecting more precise measurements 
of the Higgs data. The distributions of chi-square in the plane of {AS"', AS^) and the 
corresponding (C^,Cg) plane are shown in Fig. [T]^a) and (b), respectively. The quantities 
C-y and Cg are very close to physical observables, and so their best values are 

a, - 1.14, Cg - 0.93 . 

The chi-square per dof for this two-parameter fit is 0.88 compared to 0.86 of the SM. The 
p-values are very similar. 

Again, we further allow AFtot to vary simultaneously with AS"'' and AS^, the best fit is 

AS'' = -0.96+°!^, A5« = -0.040l°:i2g, AFtot = 0.027tlil MeV, xVdof = 17.55/19 = 0.92. 

The 2-dim contours for the correlations among the 3 parameters are shown in Fig. |2} Note 
that anticorrelation between and Cg shown in Fig. ^h) is modified to the shape shown 
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in Fig. [2](d). The elongation along the Cg is allowed with the increase in AFtot such that 
the production in ggF increases but the decays in various channels decrease. It is clear that 
including AFtot in the fit does not improve the fit. 

It is easy to notice that AS"' and AS^ were very efficient in reducing the of Higgs 
data before the Moriond update, because both the ATLAS and CMS had the diphoton data 
on the excess side of the SM value. However, after the Moriond update the CMS diphoton 
data is smaller than the SM value while the ATLAS is still larger, and therefore the 
cannot be reduced effectively no matter how AS"'' and AS*^ are varied. 

D. Vary C^, C^, Cf , C^, while keeping AS^ = AS^ = 

This choice is motivated in the scenario where there are no new particles running in the 
triangular loops of Higgs boson decaying into gg or 77, but only modifications of Yukawa 
couplings. It can be realized in a two-Higgs doublet model with a very heavy charged Higgs 
boson and modifications of Yukawa couplings can be expressed in terms of the mixing angle 
a and tan (3. 

First, we notice that there is an overall symmetry: 

simply obtained by flipping the overall sign in Eqs. ([T]) and ([2]). Furthermore, from Eq. ([s]) 
the diphoton production rate depends on \S'^\^ + and so only the relative signs of 

gauge and Yukawa couplings are important. Therefore, in the following we fix the sign of 

to be positive, while the other 3 parameters can be either negative or positive. 

Since the contributions of the bottom and the charged-lepton sectors to the diphotons are 
very small, we expect an approximate symmetry: O and Cf O . Even if we 

change the best-fit point by Cf — t- — Cf the total only changes by O(+0.01). This simply 
means that changing the sign of bottom- and charged-lepton Yukawa couplings would not 
affect significantly the loop contributions of iJ77 and Hgg. We show in Fig. [3] the 2-dim 
contours for the correlations of any 2 of the 4 parameters (Cf , , Cf , C^,). In the figure, 
we can see an approximate symmetry: o ^C'f ^-nd Cf o —Cf. Note that this figure is 
for the data after the Moriond update. 

On the other hand, the sign of is important, as shown by the two islands in the panel 
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(a). It is well known that the SM W boson and top quark loop contributions to H'j'j come 
in opposite sign. Therefore, by flipping the sign of the top quark contribution (Cf — j- — ) 
it can enhance the H'yy vertex. 



1. Before Moriond 



Since before the Moriond update both the CMS and ATLAS diphoton data are in excess, 
the dynamics of the fit indeed prefers < for positive Cy, shown in the fifth column of 



the upper half of Table VI In this way, the diphoton rate is pushed up to fit well with the 
data and significantly reduces the x^- Thus results in x^/dof = 10.46/18. 



2. After Moriond 

Nevertheless, the new CMS diphoton data affect the fit significantly. The dynamics of the 
fit cannot force the parameters to go into one direction to reduce the x^, because the ATLAS 
diphoton data is on the opposite side of the CMS data. Thus, the top- Yukawa ~ 0.8 



(see the fifth column of the lower half of Table VI), which means that the top contribution 
to the H'j'y vertex is only reduced by a small amount. Therefore, we only obtain an overall 
x'^/dof = 17.82/18 = 0.99, which is worse than the SM fit. 

In Fig. |4f^a), we show the corresponding confidence-level regions in the {C^,Cg) plane. 
The central values are Cy = 1.09, Cg = 0.91. Note that original anticorrelation between Cy 
and Cg shown in Fig. is modified to the shape shown in Fig. ^a). The enlargement 
region in Cg can be understood when C^ increases, Cg will increase but Cy is reduced such 
that the diphoton rate stays about the same. Another enlargement region in C^ direction 
can be understood as the left island of Fig. [sj^a), in which C^ is negative, such that C^ is 
large and Cg is about the same. In part (b), we show the correlation between C^ and Czy 
There are 2 islands corresponding to those shown in Fig. ^a). Both Cj and Cz-y increase or 
decrease in the same direction, though the enhancement in Czy is always smaller than Cy. 
At the best-fit point, we find Cz-y = 1-05. 
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E. Vary C^, C|, , C,, AS^ 



In this fit, we group these 6 parameters into 2 sets: (Cf , Cf, Cf , C^) and (AS"*", AS^). 
We first show the correlations among the first set in Fig. |5} which can be compared directly 
to the corresponding panels in Fig. [3j It is easy to see that all the confidence-level regions 
are enlarged due to more dof in (AS^, AS^). The next correlations between (Cf , C„) and 
(AS"'',AS'^) are shown in Fig. [oj The correlations between (Cf , Cf) and (AS"^,AS'^) are 
shown in Fig. [7j The corresponding confidence-level regions in (A5"'', AS^) and (C-y, Cg) 
planes are shown in Fig. [8| Note that these figures used the data after the Moriond update. 



1. Before Moriond 



The fit before Moriond is shown in the last column of the upper half in Table |VI[ The 
most significant changes from the previous fit are the widening of and the shift of the 
best value of from approximately —0.9 — ?■ 0. We shall explain it shortly. The and 
remains approximately the same. Although the Cf flips the sign, the absolute value 
is about the same. The flipping of the sign of Cf is simply a numerical artifact that the 
difference in is only O(10~^). In this 6-parameter fit, the signs of the best values of C^, 
Cf , Cf , are all positive, in accord with the SM. 

The shift of from —0.9 — )■ can be understood from the numerical expression of S"' 
in Eq. ([s]). In the 4-parameter fit where A5"^ = 0, the top- Yukawa coupling is made 
negative in order to increase numerically the S"^; whereas in the 6-parameter fit the goes 
to zero and A5"'' goes to a negative value to enhance S"' . This explains the shift of and 
anti-correlation between and AS'^' . 

The resulting 'X^ /dof = 9.89/16, which is pretty good. The dynamics of the fit raises the 
diphoton rate to fit the data well, so that the is reduced substantially. 



2. After Moriond 



As shown in the last column of the lower half in Table VI, the changes from approx- 
imately 0.8 — > 0, while the Cf , , and remains about 1. The Yukawa couplings are 
almost in good accord with the SM, except for the top- Yukawa. Instead, AS"'' and AS^ 
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become nonzero to accommodate the data. Nevertheless, the reduction of is very smalL 
The resulting 'X^ /dof = 16.89/16, which is worse than all other fits. 

The next correlations between (Cf , C^) and {AS"' ,AS^) are shown in Fig. |6} Both the 
AS^ and AS^ are anti-correlated with C^. On the other hand, AS"'' and AS^ do not 
correlate with C^. Similarly, the correlations between (Cf, C/) and (AS"'',AS'^) are shown 
in Fig. [7| which are also negligibly correlated. 

The corresponding confidence-level regions in (AS"^, AS^) and (C^, Cg) planes are shown 
in Fig. |8} The prediction for Cz-y is also shown. Again, the Cz-y increases or decreases in 
the same direction as C^, but is always smaller than C^. 



F. Concluding remarks 

The best-fit values for various CP-conserving fits using the Higgs data before the Moriond 



2013 are shown in the upper half of Table VI while using the data after the Moriond are 



shown in the lower half. We also show the p- value of each fit. It is clear that most of the 
fits have better p-values than the SM one before the Moriond; while all the fits are worse 
than the SM one after the Moriond. 

Before the Moriond update, the diphoton signal strength 77 dominates 

the chi-square. Both the ATLAS and CMS diphoton data are on the same side of excess 
of the SM value. The signal strength of — t- — ?■ 77 depends largely on S"^ and 5*^, 
which in turns depend mostly on and AS"^. Indeed, we have shown in the 4-parameter 
analysis and in the 6-parameter analysis, the is mostly sensitive to and AS"''. In the 
4-parameter analysis where AS"'' = 0, the best-fit value of ~ —0.9 in order to enhance 
S"'; whereas in the 6-parameter analysis ~ and AS''' ^ — 1.2 are preferred in order to 
enhance S"'. The value for AS^ also changes according to the change in the value of as 
in the expression for in Eq. ([s]). 

The and AS"'' are the two parameters most sensitive to the signal strength of pp — )■ 
H — )• 77. Keeping the other parameters as = Cf = Cf = 1, AS*^ = 0, AFtot = 0, we 
vary and AS''' only and we find the best-fit values are 

C^ = 0.92t'of,l, AS^ = -2.62tlZ xV do f = 11.17/20 . (36) 
This is the best y^/dof that we found when only two parameters are allowed to vary, 
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although the total is only 0.1 unit better than the fit with AS"^ and AS*^, which is 
statistically insignificant. 

With the Higgs updates during the Moriond 2013 meetings |24j the uncertainties in most 
channels are reduced. The decay channels other than the diphoton also began to play 
important roles in the global fits. The most dramatic change is the CMS diphoton data, 
in which the central value (the untagged) changes from 1.42 to 0.78. Now the CMS and 
ATLAS diphoton data are on the opposite side of the SM value. The dynamics of the fit 
cannot do anything to effectively reduce the from the diphoton data. We found that all 
the fits give a p-value worse than the SM one. 



V. CP VIOLATING FITS 

We devote this section to including the pseudoscalar Yukawa couplings and the pseu- 
doscalar contributions AP'^ and AP^. 



A. AST, ASS, AP^' and AP^ 
1. Before Moriond 

We have learned from all CP-conserving fits in the last section that the most efficient 
parameters fitting the data are the deviation AS"^ to the H'-fy vertex and the up-type 
Yukawa coupling , as well as the corresponding deviation AS^ to the Hgg vertex. In 
order to understand the effects of pseudoscalar nature of the Higgs boson, we first perform 
the analysis by varying the scalar contributions AS"^ and AS^, as well as the pseudoscalar 
contributions AP'^ and AP^ to the H'j'j and Hgg vertices. We keep all other parameters 
at the SM values, = C| = Cf = ^ = 1, = = C[ = and AFtot = 0. 

The best-fit parameters and the corresponding for this case are shown in the second and 



third columns of the upper half in Table VII The total = 11.26, almost the same as the 
total = 11.27 of the case varying AS"', AS^ only. Therefore, including the pseudoscalar 
contributions does not improve the fit at all. In fact, the x^ /dof is worsened. 
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2. After Moriond 



The confidence-level regions in the [AS^ , AP^), in the (A5"'', AP'''), and in the corre- 
sponding {C^,Cg) planes are shown in Fig. [9] The nearly-physical values for Cy ~ 1.1 and 
Cg ~ 0.9, which are the same as the fit using just and AS^. In order to understand 
the behavior shown in Fig. [9]^a) and (b), we can use the numerical expressions for S^, P"' in 
Eq. (g and S'^,P'^ in Eq. g. 



Numerically, 



1.1 



;-6.64 + A5^)2 + (APT)2 
(-6.64)2 



Cr, ^ 0.9 



(0.65 + A5f )2 + (APa) 



(0.65)2 

Therefore, we obtain 2 ellipses 

(7.3)2 _ g4 ^ ^^7)2 ^ (^p7)2 ^ 

(0.59)2 _ ^Q_g5 ^ ^^9^2 ^ (^p<;)2 _ (37) 

that explain the ellipses shown in Fig. [9]^a) and (b). It is clear that nonzero values of AP''' 
and AP^ are not ruled out, and the data allow for both scalar and pseudoscalar values in 
special combinations for H'j'y and Hgg vertices. 

The best-fit parameters and the corresponding for this case are shown in the second 
and third columns of the lower half in Table VII The total = 17.55, the same as the case 
varying AS"'', AS^ only. Again, including the pseudoscalar contributions does not improve 
the fit at all. 



B. Ci, CP and 

From the results of the last section, we can see that among all the Yukawa couplings 
the fitting result is more sensitive to the up-type Yukawa couplings C^. This is easy to 
understand because the diphoton data currently dominate the chi-square, and the top- 
Yukawa and HWW couplings are the most important to determine the diphoton signal 
strength. Without loss of generality we only consider the up-type scalar and pseudoscalar 
Yukawa couplings and C„ in this subsection. The effects of down-type and charged-lepton 
pseudoscalar couplings are similar but much milder. 
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1. Before Moriond 



The best-fit parameters and the corresponding for this case are shown in the fourth 
and the fifth columns of the upper half in Table VII The total = 10.53, almost the 
same as the total x^ = 10.46 of the case varying C^^^ and C^. Therefore, including the 
pseudoscalar Yukawa coupling does not improve the fit. Nevertheless, the p-values are 
about the same, and so the data cannot rule out the combination of scalar and pseudoscalar 
Yukawa couplings. 



2. After Moriond 

The 2-dim confidence-level regions among the parameters {C^,C^ ,0^) are shown in 
Fig. 10 We can directly compare Fig. [sj^a) and Fig. 10 a). The 2 islands in Fig. [sj^a) 



are now linked together in Fig. |lO[a), due to the variation of an additional parameter C^. 
The sickle-shaped region in part (c) indicates that and satisfy some equations of 
ellipses. From the numerical expressions for S^, in Eq. Q and S^, pa in Eq. (g), we have 

(7.3)2 ^ (_g 4 ^ ^ 7g ^5)2 ^ (2.78 C^f 

(0.59)2 _ (0.688 Cf)2 + (1.047 Cf)2 (38) 



which can then explain the shape in part (c). 

The correlation between and Cz-^ is shown in part (e) and that for the CP-violating 
observables, which are proportional to 2CfC^/(Cf^ + C^^), is shown in part (f). The 
Cz-y increases and decreases in the same direction as but always smaller than C^. At 
the best-fit point, ~ 1.1 while Cz-y ~ 1.05. The CP-violating observables arised from 
the mixing between scalar and pseudoscalar contributions are in general proportional to 

The best-fit parameters and the corresponding x^ for this case are shown in the fourth 
and the fifth columns of the lower half in Table |VII The p- value is slightly better than the 
CP-conserving case of varying ^ and 
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VI. DISCUSSION 



In this work, we have estabhshed a model-independent framework that enables one to 
analyze all the observed Higgs boson signal strengths and to fit to the Higgs boson couplings 
to fermions, W/Z bosons, 77, and gg. In the future when the Z7 data are available, the 
current framework can also cover that. Right now we give predictions for Z7 signal strengths. 
We have performed global fits to all the Higgs signal strengths recorded by ATLAS, CMS, 
and at the Tevatron. Furthermore, we have also performed fits including the pseudoscalar 
up-type Yukawa coupling and the pseudoscalar contributions P'^ and in the H'yy and 
Hgg vertices, respectively. 

Furthermore, we have performed fits with respect to all the Higgs data collected before 
and after the Moriond 2013 meetings. The main reason why we separately performed them 
is because of the dramatic change in the fits due to the shift of the CMS diphoton data from 
1.42 to 0.78 of the SM value. Before the Moriond update the dynamics of the fit pushes the 
parameters to increase the diphoton rate such that the is reduced effectively. However, 
with all the data after the Moriond the dynamics of the fit cannot find the optimal set of 
parameters so that the resulting x^/rfo/'s are indeed worse than the SM. 

In summary, the p- value of the SM Higgs boson is 0.65 performed with all the data after 
the Moriond. Its p-value is higher than any other fits considered in this work, both the 
CP conserving ones and the CP-violating ones. We also plot the p-values for all the fits 



considered in this work in Fig. 11 



Our findings are summarized as follows. 

1. Before the Moriond, the SM already enjoyed a good x^/dof = 0.8 (p-value = 0.74), 
which means that the SM Higgs description of the data is reasonably well. Out of 
the total = 17.5 about one half comes from the — )■ 77 data. Similarly, after the 
Moriond the SM has x'^/dof = 0.86 (p- value = 0.65). 

2. If only the total Higgs boson width is allowed to vary, we are able to constrain the 
deviation AFtot to be less than 1.2 MeV at 95%CL. Given the SM Higgs boson width 
is about 4.1 — 4.2 MeV for Mh = 125 — 126 GeV, the nonstandard decay branching 
ratio of the Higgs boson is less than 22% at 95% CL. This is a real improvement from 
previous estimates of about 40% [82]. 
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3. The most efficient set of parameters to fit to the data before the Moriond are the 
additional particle contributions to the loop functions of if 77 and Hgg vertices, AS"'' 
and AS^ respectively. This is because both the ATLAS and CMS have diphoton 
data above the SM value. The best x^/dof obtained is about 0.56. This is easy 
to understand, as the total is currently dominated by ii — )■ 77 signal strength. 
Nevertheless, with the Moriond update the CMS diphoton data is now below the SM 
value. No optimal set of parameters can be found to effectively reduce the total x^- 

4. With the data before the Moriond, another efficient set of parameters are and AS"''. 
We found that they are equally effective as (AS'^, AS^). Effectively, the modification 
in takes up the place of AS^. Again, the reason is the domination of if — 77 in 
the total x^- Nevertheless, after the Moriond update no optimal set of parameters can 
be found. 

5. The relative HVV coupling now stands at = 1.01 IIqH in the 6-parameter fit (and 
a similar value in the 4-parameter fit). This implies that the observed Higgs boson 
accounts for most of the EWSB, and leaves little rooms for additional Higgs bosons 
that are also responsible for EWSB. Nevertheless, if we take —2a to the central value, 
the can be as low as 0.7. The vector-boson scattering could become strong if the 
UV part of the Higgs sector is very heavy [83] . 

6. The current data do not rule out the pseudoscalar contributions to the ii77 and Hgg 
vertices nor the pseudoscalar Yukawa couplings. Nevertheless, including pseudoscalar 
contributions to ii77 and Hgg vertices or pseudoscalar Yukawa couplings do not im- 
prove the fits. The current Higgs observables are not sensitive to CP-violating effects, 
and so only combinations of scalar and pseudoscalar contributions are constrained, as 



shown in Eqs. (37) and (38). Thus, the current Higgs data do not rule out or favor 



pseudoscalar couplings. 

Era of Higgs-boson precision studies now begins - Higgcision. In this work, we have 
already seen dramatic changes in the fits with the data collected before and after the Moriond 
2013. We are looking forward to more and more data in the upcoming Summer 2013 and 
the following years. 
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TABLE I. Data on signal strenghts of — t- 77 recorded by ATLAS and CMS, and at the Tevatron 
before and after Moriond 2013. The luminosity updates at 8 TeV are shown in the parenthesis. 
The percentages of each production mode in each data are given (details are given in the text). 
The of each data with respect to the SM is shown in the last two columns for before and after 
Moriond. The sub-totals of this decay mode are shown at the end. 

Channel Signal strength /j, M//(GeV) Production mode Xsm(^^^^) 



Before After ggF VBF VH ttH Before After 



ATLAS (4.8fb-i at 7TeV + 13.0 (20.7) fb"! at 8TeV): [9l[25] 


fJ'ggH+ttH 


1.8 ±0.49 


1.6 ±0.4 126.8 100% 


2.67 


2.25 




2.0 ±1.4 


1.7 ±0.9 126.8 - 100% 


0.53 


0.60 




1.9 ±2.6 


1.8t\i 126.8 - - 100% 


0.12 


0.38 


CMS (5.1fb-^ at 7TeV + 5.3 (19.6)fb-^ at 8TeV) [I31E7] 


untagged 


-1 49+0.55 


0.78l[jJ^ 125 87.5% 7.1% 4.9% 0.5% 


0.73 


0.62 


VBF tagged 


9 9c; + 1.34 


2.25+\-^l 125.8 17% 83% 


1.44 


1.44 


Tevatron (lO.Ofb^^ at 1.96TeV): ^ 


Combined 




6.14^3-25 125 78% 5% 17% 


2.60 


2.60 






subtot: 


8.09 


7.89 
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TABLE II. The same as Table I but for H ZZ^*\ 



Channel 


Signal strength /i 


Mj^(GeV) Production mode 


XsmO 


sach) 




Before After 


ggF VBF VH ttH 


Before 


After 




ATLAS (4.8fb" 


at 7TeV + 13 (20.7)^"^ at 8TeV) [TOl [26] 






Inclusive 


1.0 ±0.4 1.5 ±0.4 


125.5 87.5% 7.1% 4.9% 0.5% 


0.0 


1.56 




CMS (5.1fb-i 


at 7TeV + 12.2 (19.6) fb"^ at 8TeV) [HEEi 






Inclusive 


0.80l°:i o.9it°io 


125.8 87.5% 7.1% 4.9% 0.5% 


0.33 


0.09 



subtot: 0.33 1.65 



TABLE HI. The same as Table E but for H WW^*\ 



Channel 


Signal strength ^ M//(GeV) Production mode 


Xsm( 


each) 




Before 


After ggF VBF VH ttH 


Before 


After 


ATLAS (4.8fb-^ at 7TeV ± 13 (20.7) fb"^ at 8TeV) [13 EF] 


Inclusive 


1.5 ±0.6 


1.0 ±0.3 125.5 87.5% 7.1% 4.9% 0.5% 


0.69 


0.00 


CMS (up to 4.9 fb~i at 7TeV ± 12.1 (19.5) fb^^ at 8TeV) [l3l|29] 


0/1 jet 


77+0-27 

^- ' ' -0.25 


0.76 ±0.21 125 97% 3% 


0.73 


1.31 


VBF tag 


-0 05+°-^^ 


-0.05lj]:^5 125.8 17% 83% 


2.01 


2.01 


VH tag 


— n "^1 +2-22 
'-'•'3 ■'--1.94 


-0.3ll?J^ 125.8 - - 100% 


0.35 


0.35 


Tevatron (lO.Ofb"^ at 1.96TeV): [7] 


Combined 


0.85t°|? 


Q.Shtofi 125 78% 5% 17% 


0.03 


0.03 






subtot: 


3.81 


3.70 
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TABLE IV. The same as Table |l] but for H — )• bb. There are no updates for this channeL 



V_.- XJ.<ApXJ.XXV_>± 




A/ftj((~reVl Prorlnrtion morle 
ggF VBF VH 


ttH 






ATLAS f4 8fb"-^ 


at 7TeV + 13 Ofb"^ at STeV) l26l 






VH tag 


—0.4 ± 1.0 


125 5 - - 100% 




1.96 




CMS fuD to 5 Ofb" 


^ at 7TeV + 12.1fb""'^ at 8TeV) [13 






VH tag 


-1 q-i+0.65 
-'-•'-'J^-0.60 


125.8 - - 100% 




0.27 


ttH tag 


-0.80i?i° 


125.8 


100% 


0.73 




Tevatron 


(lO.Ofb"^ at 1.96TeV): [8] 






VH tag 




125 - - 100% 




0.59 



subtot: 3.55 



TABLE V. The same as Table |T] but for H tt. The correlation for the rr data of ATLAS is 
p = —0.50 and —0.49 before and after Moriond, respectively. The percentages of the production 
modes differ very tiny before and after Moriond. 



Channel 


Signal strength fj, M//(GeV) Production mode 


xiM(each) 




Before 


After ggF VBF VH ttH 


Before 


After 


ATLAS (4.6fb"^ at 7TeV + 13.0fb"^ at 8TeV) HUES] 




2.38 ± 1.57 


2.30 ± 1.60 125.5 100% 


1.60 


1.41 


niVBF + VH) 


-0.25 ± 1.02 


-0.22 ± 1.06 125.5 - 59.4% 40.6% 






CMS (up to 4.9fb-i at 7TeV + 12.1 (19.4) ib'^ at 8TeV) [13 [30] 


0/1 jet 


n oc+0.68 


0-76io|2 125 77.8% 13.8% 7.6% 0.8% 


0.05 


0.23 


VBF tag 


0-82tr75 


1.40+°|? 125 20.9% 79.1% - 


0.05 


0.49 


VH tag 




0.771};^^ 125 - - 100% 


0.005 


0.02 






subtot: 


1.70 


2.15 
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TABLE VI. The best fitted values and the la errors for the parameters in various CP conserving 
fits and the corresponding chi-squarc per degree of freedom and the p-value before and after the 
Moriond 2013. The p-values for the SM fit are 0.74 and 0.65 for the data before and after the 

Morioud. respect iv(^ly. 





Vary APtot 


Vary AS^, 


Vary AS^, 


Vary C|, 


Vary C^, C^, C^, Cy 


Parameters 




ASa 


AS9, APtot 


Cf , Cy 


AST, AS9 


Before Moriond 




1 


1 


1 


-0.88l°i? 


0.00 ± 1.13 


cl 


1 


1 


1 


1 -19+0.45 


1.19l°il 


cf 


1 


1 


1 


-0.97l°:i 


0.98 ±0.30 




1 


1 


1 


0.97lgii 


^•^0-0.15 


AST 





r, 70+1.11 







1 90+2.44 
— i.^o_2.49 


ASa 





-0 n50+°°6^ 

U.UOU_Q 065 







n 70+O.8I 
'-'•'''-0.80 


APtot (MeV) 


-0.022l°i3 





0.79t?:?} 








xVdof 


17.48/21 


11.27/20 


10.83/19 


10.46/18 


9.89/16 


p- value 


0.68 


0.94 


0.93 


0.92 


0.87 


After Moriond 




1 


1 


1 


0.80l°;J^ 


0.00 ± 1.18 




1 


1 


1 


-0.98t°:ll 






1 


1 


1 


0.98lg;il 


1.01 ±0.23 




1 


1 


1 




i.oil°:ll 


AST 





-0.96l°:i 







0.78l|i 


ASa 





-0.043 ± 0.052 


-o.o4ot[J:i2g 







APtot (MeV) 


^•^^'-'-0.41 





0.027tlil 








xVdof 


18.89/21 


17.55/20 


17.55/19 


17.82/18 


16.89/16 


value 


0.59 


0.62 


0.55 


0.48 


0.39 
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TABLE VII. The best fitted values and the la errors for the parameters in the CP-violating fits 



and the corresponding chi-square before and after Moriond 2013. 



Parameters 


Vary AS^, A5»,APT, AP^ 


Vary C^C^ , 


Before Moriond 




1 


1 


r\ r-^+0.65 


n t; ,1+0.65 
-0-54Io.39 


Cd 


1 


1 


1 


1 


cf 


1 


1 


1 


1 




i 


i 


n Qo+o.io 

U.»0_Q^2 


n Qo+o.io 

U.»0_o 12 


AS'''' 


r> 00+16 9 


n OQ+17 4 
-U.^»_3g0 








ASS 


f;9+0-63 
^•0^-0.70 


— n «9+0-63 








APtot (MeV) 






















-n 46+^-^^ 

"■^'J-0.22 


'-'•^'J-0.22 






-6.3li^«,« 










o.6ot?:°f 


o.6ol?:°f 








xVdof 


11.26/18 


11.26/18 


10.53/19 


10.53/19 


value 


0.88 


0.88 


0.94 


0.94 


After Moriond 




1 


1 


0.48t°:|^ 


0.48l°i^ 


CI 


1 


1 


1 


1 


cf 


1 


1 


1 


1 




1 


1 


0.995lHgI 


0.995l°:?gl 


AST 


-0.921^^3°° 


0.7ll^^5f 








AS9 


-0.55^°:^^ 


-0.6<^^ 








APtot (MeV) 






















U.0U_Q4Q 


-o.5ol°:tt 


APT 


n 77+7-67 


4 7c;+3.69 
^•''3_i3.i9 








APs 


^■"^-0.06 


'-'•"-'--0.052 








xVdof 


17.55/18 


17.55/18 


17.17/19 


17.17/19 


]9- value 


0.49 


0.49 


0.58 


0.58 
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FIG. 1. The confidence- level regions of the fit by varying AS"*" and AS^ only, (a) in the {AS^ , AS^) 
plane and (b) in the corresponding (C^, Cg) plane. The contour regions shown are for Ax^ < 2.3 
(red), 5.99 (green), and 11.83 (blue) above the minimum, which correspond to confidence levels of 
68.3%, 95%, and 99.7%, respectively. The best-fit point is denoted by the triangle. 
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FIG. 2. The confidence-level regions of the fit by varying AS"'', A5^, and AFtot, (a) in the 
(A5"'', AFtot) plane, (b) in the (AS^, AFtot) plane, (c) in the (A5'^, AS^) plane, (d) in the corre- 
sponding (C^y, Cg) plane. The description of contour regions is the same as Fig. [TJ 
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FIG. 4. (a) Same as Fig. [s] but in the corresponding plane of {C^,Cg). (b) Prediction in the 
corresponding (C^, Cz-y) plane. The description of contour regions is the same as Fig. [Tj 
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FIG. 5. The confidence-level regions of the fit by varying C^, Cf , Cf , C^,, AS"'' and AS^ while 
keeping AFtot = 0. Shown are the correlations among (C^, C^, , Cy). The description of 
contour regions is the same as Fig. [!} 
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FIG. 7. Same as Fig. [5] Shown are the correlations between (C^, Cf) and (AS"^, A5^). 
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FIG. 8. Same as Fig. [5] Shown are the correlations between I^S'^ and A5^, between and Cg, 
and between and 
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FIG. 9. The confidence-level regions of the fit by varying the scalar contributions AS"' and 
AS^, and the pseudoscalar contributions AP"/ and AP^ while keeping = = Cf = 1, 
= = Cf = and AFtot = 0. The description of contour regions is the same as Fig. fl| 
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FIG. 10. The confidence-level regions of the fit by varying the scalar Yukawa couplings and 
C^,, and the pseudoscalar Yukawa couplings C^; while keeping C'^ = Cf = 1, = Cf = 0, and 
AS"'' = = /S.P'^ = AP^ = AFtot = 0. The description of contour regions is the same as Fig. [T| 
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No. of fitting parameters 

FIG. 11. The p-values for various fits considered in this work, including CP-conserving (CPC) 
and CP-violating (CPV) ones. The CPC cases with the No. of fitting parameters equals 1 denotes 
the case varying only AFtot! 2 denotes varying only AS'^ and AS^; 3 denotes varying only AS'''', 
AS^, and AFtot; 4 denotes varying only C^, i and Cy] 6 denotes varying C^, C^, Cy, 

AS"', and A5^. The CPV cases with the No. of fitting parameters equals 3 denotes varying only 
C^, C7^, and C^; 4 denotes varying only A5^, A^^', AP^, and APf . 
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